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Abstract 

Using multicanonical simulations, we study the effect of charged end groups on helix formation 
in alanine based polypeptides. We confirm earlier reports that neutral poly-alanine exhibits a 
pronounced helix-coil transitions in gas phase simulations. Introducing a charged Lys"*" at the 
C-terminal stabilizes the a helix and leads to a higher transition temperature. On the other hand, 
adding the Lys"*" at the N-terminal inhibits helix formation. Instead, a more globular structure was 
found. These results are in agreement with recent experiments on alanine based polypeptides in gas 
phase. They indicate that present force fields describe accurately the intramolecular interactions 
in proteins. 
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I. INTRODUCTION 



First principle calculations of protein folding have remained a computationally hard prob- 
lem. This is in part because the energetic and entropic balance of a solvated protein in 
equilibrium is complicated to express on a computer. Other problems arise from the in- 
herent difficulties in sampling the rough energy landscape of a protein. It is not a priori 
clear what contributes most to the numerical shortcomings: poor sampling techniques, lack 
of accuracy in the force fields that describe the intramolecular forces in a protein, or the 
crude approximations in the modeling of protein-solvent interaction. The introduction of 
high-resolution ion mobility measurements, which allow experiments to be performed on 
biological molecules in gas phase , has opened a way to simplify the problem. By com- 
paring simulations with gas-phase experiments one can avoid the complications arising from 
the modeling of protein-solvent interactions. Hence, gas-phase experiments offer a way to 
test whether current force fields describe correctly the intramolecular forces and the ability 
of present simulation techniques to sample accurately low-energy configurations. 

Such evaluation of methods and energy functions through comparison with gas phase 
experiments is done best with simple systems. Often used are homopolymers of amino 
acids. For instance, the helix-coil transition has been studied extensively with poly-alanine 
0, y. Is]. Numerical results indicate a strong propensity to form helices if the temperature 
is below a certain critical value. On the other hand, the experimental observations seemed 
to contradict those theoretical results. For this reason, it was conjectured in Ref. gI that 
the strong helix-coil transition for polyalanine in gas-phase simulations is an artifact of 
the utilized energy function. However, the deviations may also result from differences in 
the systems studied. Gas phase experiments require partially charged molecules while in 
simulations it is commonly assumed that the molecule has no total charge. 

In order to settle this question and quantify the effects of charged end groups we have 
simulated three alanine based polypeptides ingas phase. The first one is Alain which we 

n p n II 

are familiar with from earlier work 17|, 18|, 19|. Simulations of this neutral molecule are 
supplemented by those of two charged polypeptides. In Ace-Lys+-Alaio the charged group 
sits at the N terminal while in Ace-Alaio-Lys^ the charge sits at the C-terminal. Our results 
confirm the experimental findings that capping the N terminal by a positively charged Lysine 
destabilizes the helix and leads to a more globular low-temperature structure, while a C- 
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terminal Lys"*" can stabilize the a helix. 



II. METHODS 



Our simulations utilize the ECE 



iPP/3 force field tldl as implemented in the 2005 version 
of the program package SMMP [UJ, [12] . Here the interactions between the atoms within the 
homopolymer chain are approximated by a sum Eecepp/z consisting of electrostatic energy 
Ec, a Lennard- Jones term Elj, hydrogen-bonding term Ehb and a torsion energy Etot'- 

-E-ECEPP/s = Ec + Elj + Ehb + Etot 
_ ^ 332gigj 

_l_ [ ^ji _ ^ij \ 

I „12 ^10 ) 

+ J2Ui{l±cos{ni^i)) , (1) 
I 

where Vij is the distance between the atoms i and j, C,i is the l-th torsion angle, and energies 
are measured in Kcal/mol. 

The above defined energy function leads to a landscape that is characterized by a multi- 
tude of minima separated by high barriers. In order to enhance sampling we therefore utilize 
the multicanonical approach 13|, as described in Ref. [l5|. Configurations are weighted 
with a term wmu (.) detenn,,;:; ite^tivelyQ such tKat the p.^abil.t, d.tnbut.o. obey. 

Pmu{E) oc n{E)wMu{E) ^ const , (2) 

where n{E) is the spectral density of the system. Thermodynamic averages of an observable 
< O > at temperature T are obtained by re- weighting fl7 |: 

J dx 0{x)e-^(^y'-^ /wMujEjx)] 
Jdx e-^(-)ABr/^^^ ^E{x)] 

where x counts the configurations of the system. 

Various quantities are measured during simulations for further analysis. These include 
the energy E, the radius of gyration Rgy and the end-to-end distance Ree as measures of the 
geometrical size, and the helicity. Here we define the helicity by the number of backbone 



hydrogen bonds, Ub-HB, between residues + 4), which is characteristic for a hehces. 
We also monitored the number of non-backbone hydrogen bonds, Unb-HB, as a measure of 
tertiary interactions. Finally we recorded the lowest energy configurations. 



III. RESULTS AND DISCUSSIONS 

As we are interested in the helix-coil transition in alanine-based polypeptides, we start 
our investigation with an analysis of the average helicity for the three peptides, measured by 
the number of backbone hydrogen bonds (rib^HB)- This quantity is displayed as a function 
of temperature in Fig. [H For the neutral molecule, Alaio, a sharp transition is observed 
that separates a disordered high temperature phase with vanishing helicity from an ordered 
low-temperature phase where most configurations have high helicity. Not surprisingly, the 
minimal energy configuration exhibits a fully formed a helix and is displayed in Fig. UK- 
Displaying the specific heat as a function of temperature in Fig. [2], we observe a sharp peak 
which indicates that this helix-coil transition takes place at a "critical" temperature Tc = 
462K ± 5K . Note that this value is about 40 K higher than the corresponding temperature 
427K ± 7K for the older ECEPP/2 force field [ol although the two transitions differ little 
otherwise (data not shown). 

In order to study the effect of charges on the stability of an a helix, we perform a second 
simulation with Ace-Alaio-Lys"*" . In this molecule, a positively charged Lys residue is added 
at the C-terminal end of the chain. The measured values for the helicity of Ace-Alaio-Lys"*" 
are similar to Alaio, forming again a single extended helix upon lowering the temperature 
(see Fig. [I]). Note that the peak in specific heat (Fig. [2]) is with Tc= 557 ± 19K at an ^ 100 
K higher temperature, and is much wider. This increase in the transition temperature 
indicates that the charged Lys"*" at the C-terminal favors the a helix through stabilizing 
the helix dipole. Parts of the a helix survive even up to higher temperatures widening the 
transition regime considerably. The stable a helix structure is also seen in the minimal 
energy configuration displayed in Fig. Hh. Note the additional hydrogen bond formed by the 
charged end group. Our results are consistent with recent experiments that also indicate 
stabilization by charged groups capping the C terminus and by the interaction of the charge 



with the helix dipole 



18| 



The third peptide that we have studied is Ace-Lys"'"-Alaio. Here the charged Lys is added 
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at the N-terminal. For this peptide one finds only a small peak in the specific heat (see 
Fig. [2]), just below that of Alaio (^c= 449 ± 37K), and a shoulder at higher temperatures. 
At low temperatures, its helicity (displayed in Fig. [T]) is substantially lower than for the 
previous two peptides. Fig. SJ: displays the minimal energy configuration. Only a small 
partial helix is seen in the central region, while the end groups coil and form a more compact 
globular structure. The inlay in Fig. [1] indicates that non-helical hydrogen bonds, stabilizing 
the globular structure, are formed already around 600K, above the small peak in spec. heat. 
This formation of the globular structure as well as the joining up of the end groups is 
corroborated by the collapse of the molecule in that temperature range, as seen in the radius 
of gyration, (Rgy), as well as in the end-to-end distance, (Ree), in Fig. [HJ For comparison, 
in the uncharged molecule the drop in (Rgy) occurs at lower temperatures, i.e. at the helix- 
coil transition, with practically no change in the end-to-end distance (Ace-Alaio-Lys"'", the 
charged helix-forming variant, shows a similar behavior). In contrast, the formation of the 
helical sub-segment in Ace-Lys^-Alaio does not show up significantly in the geometrical 
measures of that molecule. Hence, our data indicate that there is no well-defined helix-coil 
phase transition for Ace-Lys^-Alaig. Instead, the shoulder in the specific heat rather marks 
the collapse of the molecule, while the peak at lower temperatures marks the formation of 
the helical sub-segment. The reason for the suppressed and only partial helix-formation is 
the unfavorable interaction between the N-terminal Lys"*" and the helix dipole, destabilizing 
the extended helix. This conjecture is again supported by experimental finding in gas phase 
where the relative cross sections of this peptide also indicate configurations that are more 

n 

compact than helices |6|]. 



IV. CONCLUSION 



This study has been motivated by a discrepancy between gas-phase experiments and first 
principle calculations of poly-alanine that either could indicate shortcomings of our energy 
function, the ECEPP force field, or differences between the systems studied. As the orig- 
inal work in Refs. ^ focused on neutral molecules, while the experiments studied charged 
molecules jsl, we have now compared simulations of polyalanine with either no charge, the 
charge at the C-terminus, or at the N-terminus. Our simulation for the charged molecules 
are now in agreement with the experimental results. Our results underline two points. First, 
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when comparing gas phase experiments with simulations one has to make sure that the ex- 
perimental settings are adequately described in the computer experiment. Secondly, the 
agreement between gas phase experiment and simulation suggests that the current genera- 
tion of force fields describes the intramolecular forces within a protein sufficiently accurate 

ten observed failure of protein 



19|) therefore result likely from 



to allow for correct secondary structure formation. The of 
simulations in finding the correct structure (see, e.g. Ref. 
insufficient sampling and poor representations of the protein-solvent interaction. 
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Figure captions: 

Fig. 1: Average helicity, measured by the number of backbone hydrogen bonds, {rih-^HB), as 
function of temperature T for the three polypeptides, obtained from a multicanonical 
simulation with 5 x 10^ sweeps, (open circles) Alaio, (open squares) Ace-Alaio-Lys"^, 
and (filled squares) Ace-Lys^-Alaio; error bars are included and are mostly about the 
symbol size or less. The inlay shows the non-backbone hydrogen bonds {unb-HB), a 
measure of tertiary interactions. 

Fig. 2: Specific heat C(T) as function of temperature T for the three polypeptides, obtained 
from a multicanonical simulation with 5 x 10^ sweeps. Symbols as in Fig. 1. Error 
bars are included and are mostly about the symbol size or less. 

Fig. 3: Average radius of gyration {Rgy) as function of temperature T for the three polypep- 
tides, obtained from a multicanonical simulation with 5 x 10^ sweeps. Symbols as in 
Fig. 1. Error bars are included and are mostly about the symbol size or less. The 
inlay shows the average end-to-end distance, (Ree)- 

Fig. 4: Minimal energy configurations of (a) Alaio, (b) Ace-Alaio-Lys"'", and (c) Ace-Lys"*"- 
Alaio, obtained by minimizing the lowest energy configuration found during the sim- 
ulation. The dashed lines show hydrogen bonds. The pictures have been drawn with 

vmdQ. 
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